Viscous fingering in volatile thin films 
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A thin water film on a cleaved mica substrate undergoes a first order phase transition between two 
values of film thickness. By inducing a finite evaporation rate of the water, the interface between the 
two phases develops a fingering instability similar to that observed in the Saffman- Taylor problem. 
We draw the connection between the two problems, and construct solutions describing the dynamics 
of evaporation in this system. 



The processes of evaporation and dewetting of thin liq- 
uid films are becoming increasingly important in modern 
technologies e.g. the fabrication of electronic chips, mi- 
crofluidic devices, and biosensors. These noncquilibrium 
processes depend on various details such as the temper- 
ature distribution in the sample, the nature of the in- 
teraction between the substrate and the fluid, and the 
substrate roughness. Depending on these details, a thin 
liquid film may, for instance, break into isolated droplets, 
evaporate uniformly, or exhibit a fingering instability 
similar to that observed in the Saffman- Taylor problemfl] 
where an inviscid fluid penetrates a viscous fluid. This is- 
sue naturally generated extensive experimental and theo- 
retical studies - see Ref.0 for a recent review of the field. 
However, due to its complexity, most theoretical stud- 
ies of the problem are based on numerical computations. 
In this paper we draw the connection between the evo- 
lution of evaporating thin films and the Hele-Shaw type 
dynamics. This connection allows us to construct analyt- 
ical solutions for the evolution of evaporation. Examples 
for these solutions are depicted in Fig. 1. 

The experimental system showing this kind of 
behavior Q, consists of a clean mica substrate covered 
by a thin water film. Under proper conditions of tem- 
perature and pressure, two films of different thickness 
(approximately 2 and 12 nm) can coexist. This behavior 
is a consequence of the competing van der Waals and po- 
lar surface forces between the water and the substrate [J]. 
A first-order phase transition between the film heights is 
induced by changing the vapor pressure. On the scale 
of nanometers, gravitational effects are negligible, and if 
the evaporation rate is sufficiently small one may neglect 
changes in the film thicknesses due to pressure gradients. 
This implies that the dynamics of evaporation, in this 
system is, effectively, that of the interface between the 
phases. 

We begin by constructing the equations of motion, and 
from now on refer to the film domain of higher thickness 
as a " droplet" . Owing to the small thickness of the wa- 
ter film (in both phases), the fluid dynamics is viscous, 
and the lubrication approximation with no-slip boundary 
conditions at the substrate- liquid interface and zero shear 
at the liquid-air interface may be employed. Thus, the 
velocity of the fluid along a profile perpendicular to the 
substrate, is approximately parabolic. Averaging over 




FIG. 1: The evolution of evaporating droplets: Left-hand 
panel shows snapshots of the contour of a droplet with five- 
fold symmetry. Right-hand panel shows a dry patch pene- 
trating, from right to left, into the liquid domain. In both 
cases successive contours differ by a fixed time interval. 

this profile one obtains Darcy's law, 

h 2 

v = --VP, (1) 

of 

where h is the film height, v is its viscosity, P is the 
pressure, and v is the (two-dimensional) velocity field. 
The dynamics outside the droplet may be ignored since 
the fluid velocity in these regions is approximately zero, 
and evaporation is negligible due to the strong adhesion 
force between the substrate and the fluid. Thus evap- 
oration takes place, effectively, only within the droplet 
where it is assumed to be uniform. Assuming the film 
height to be approximately fixed, the variation in the 
two-dimensional fluid density, p, may be also neglected 
and therefore uniform evaporation implies that 

pV • v = -k, (2) 

where k is the rate of evaporation (in units of mass per 
unit time per unit area) . Combining the above equations 
one concludes that, inside the droplet domain, T>, the 
pressure satisfies Poisson's equation, while outside this 
region, it is constant. Without loss of generality we set 
this constant to be zero, and thus 

V 2 P = /3 within V, , . 

P = otherwise, 

where (5 = 2>vn/{h 2 p). The above equation describes the 
pressure away from the droplet's interface. Near that 



2 



interface, additional degrees of freedom, such as the film 
height and the boundary curvature, come into play. The 
local behavior near the droplet 's edge sets the boundary 
conditions for the above equations. In the regime where 
the surface tension energy associated with the interface 
can be neglected, the pressure is continuous and therefore 
the boundary condition for Eq. ([3]) is 



P(z) = G for z € dV, 



(4) 



where, dT> denotes the interface contour. Hereafter 
z = x + iy will represent a complex coordinate on the 
mathematical complex plane where the droplet resides. 
Eqs. ^ and (@| are equivalent to those of the "stamping 
problem" [5|, ||| . The latter is associated with the evolu- 
tion of droplets in Hele-Shaw cell when the gap between 
the two parallel plates changes in time . 

Important characteristics of the droplet 's shape are 
provided by the set of the (interior) harmonic moments: 



tk = — / dxdyz k . 
n Jzev(t) 



(5) 



where k is a non-negative integer, and the integration 
is over the droplet's domain (here, for clarity, we indi- 
cate its explicit time dependence). Following Entov and 
Etingof [5| , we take the time derivative of the above equa- 
tion, 



dt 



1 



— — = — J> dlz k Vr, 
dt 7T J z £Q V ( t ) 



(6) 



Here v n = — (h 2 /3v)d n P is the interface velocity, where 
d n is the derivative normal to the interface. Since P, 
vanishes on the boundary, one may add to the integrand 
a term of the form (h 2 /2>v)Pd n z k . This enables us to use 
Green's theorem, 

^T = ^rJ dl(Pd n z k -z k d n P) 



dt 



37W 



dV(t) 



v(t) 



dxdy (PV 2 z k - z k V 2 P) 



(7) 



However, analyticity implies that V 2 z k = 0, and substi- 
tuting Eq. ([3]) one immediately obtains dtk/dt = — ifc/V. 
Thus all tk decay at the same rate, 



tk=tk{0)e '/ T , where T = p/n. 



(8) 



Knowing the precise time dependence of the harmonic 
moments allows one to construct exact solutions for the 
evolution of the droplet. In particular, the exponential 
decrease of t , which is proportional to the droplet's area, 
implies that the interface velocity of a circular droplet is 
proportional to its radius, in agreement with the experi- 
mental data 0. 

It is instructive to consider the evolution of a droplet 
whose initial shape is an ellipse. The harmonic moments 




FIG. 2: The evolution of an evaporating elliptical droplet 



in this case vanish for all k > 2, and since Eq.© implies 
no generation of new moments, the shape of the droplet 
at any later time remains elliptical. Thus, the pres- 
sure takes the form P = (/3/4)Re [\z\ 2 - a{t)z 2 - j(t)] , 
where a(t) and j(t) are time-dependent functions to be 
determined (linear terms in x and y are omitted since 
they lead merely to a trivial shift of the droplet). This 
form of P satisfies Poisson's equation ([3]), while the 
boundary condition (j4|) is satisfied on an ellipse, 



-(*) 



-(*) 



1. 



(9) 



provided \a(t)\ < 1 and j(t) > 0. Here a± are the semi- 
major and semiminor axes of the ellipse, 



a 2 ± 



7(f) 

it <*(*)' 



(10) 



where without loss of generality we choose a(t) to be 
real and positive. From Darcy's law ([1]) one can evaluate 
the velocity of the boundary and deduce that dta± = 
— [1 T at(t)]a± / (2t) . The solution of this equation and 
(fTU| yields 



a(t) = 



1 



y/l + Ae- 2t lt 



7(0 



Be -2t/r 

y/l + Ae- 2t / =T 



(11) 



where A and B are constants determined by the initial 
shape of the droplet. Substituting pT|) into (TTU)) we ob- 
tain 



al(t) 



Be ~2t/r 



Vl + Ae- 2t / T T 1 



(12) 



The evolution of an elliptical droplet, as follows from ([9]) 
and (fl~2|) . is depicted in Fig. 2. It shows that the droplet, 
eventually, contracts to the straight line which connects 
the ellipse foci. In more general cases the droplet may 
contract to a skeleton of line segments called a " mother 
body" Q. In reality, however, surface tension effects ne- 
glected here are expected to modify the evolution at its 
late stage, and to break the narrow ellipse to a line of 
small circular droplets. 

To describe the evaporation of droplets with more com- 
plicated shapes, we shall employ conformal maps from 
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the interior of the unit circle (w domain) to the interior 
of the droplet (z domain). A simple example of such a 
map is the Joukowski map, 



z(w) = r(i)i 



v - c(t) 



•ui(t) 



(13) 



with \c(t)\ > 1. The time evolution of the droplet is de- 
termined by the time dependence of the coefficients of 
the map, r(t), c(t), Uo(t), and Ui(t) in the above exam- 
ple. This time dependence is dictated by the exponential 
decay of the harmonic moments ([5]). In order to see how, 
let us introduce the Schwarz function [ij], S(z), associated 
with the droplet's boundary in the z plane. This func- 
tion is analytic in a strip containing the boundary, and 
satisfies the condition z — S(z) for z 6 dT>. Hereafter, 
an overbar will denote complex conjugation. Expressed 
in terms of the conformal map from the interior of the 
unit circle to the droplet domain, the Schwarz function 
takes the form 



S(z) = 



1 



w(z) 



(14) 



where w(z) is the inverse map from the droplet boundary, 
in the z plane, to the unit circle in the w plane. This 
is because, on the unit circle, w = 1/w, and therefore 
S(z) = z(w) = z(w) = z(l/w). 

The relation between the Schwarz function and the 
harmonic moments of the droplet's domain is revealed 
by expressing the integral (O, using Green's theorem, as 
a contour integral along the droplet's boundary, 



t k = 



1 

2tt7 



av{t) 



dlS{z)i 



(15) 



Knowing the time dependence of the harmonic moments, 
Eqs. (Ti"4"|) and (fTS"]) , determine the time evolution of the 
conformal map z(u>). 

Consider, e.g., the Joukowski map (|13[) . where without 
loss of generality, we choose its coefficients to be real. The 
Schwarz function, in this case, takes the form 



S(z) 



r(t) u {t)w{z) 



w(z) 1 — c{t)w{z) 



ui(t). 



(16) 



Note that the inverse map iv(z), being a solution of 
quadratic equation, is defined on two Riemann sheets, 
thus S(z) is analytic on a Riemann surface of genus zero 
constructed by gluing the two patches along the cut, 
which lies outside the droplet domain. Two poles of S(z) 
are located inside T> on the physical sheet, while two 
other poles are at z — > oo of both sheets. Thus within 
the droplet domain 



3=1 



(17) 




FIG. 3: The evolution of droplet pinching due to evaporation. 
Left-hand panel shows the full droplet while the right-hand 
panel is a magnification of the area where the droplet pinches 
off. 



with no other singularities. Performing the contour in- 
tegral (fl~5"j). by the residue theorem, we obtain tk = 
Ylj^jlj- Now, since all tk decay at the same rate, q.j 
should be time independent, while Lij should decay as 
Mi(t) = ^(0)e-*/r. 

To relate the parameters of the Schwarz function, qj 
and pij, to those of the conformal map r, Uj, and c, we 
require both expressions, (fTBf and JT7J), to have the same 
poles and residues. This condition leads to the relations 



qi = z(Q), 

dz 

Mi = r — 

aw 



92 = 
M2 = 



z(l/c), 
uq dz 
c 2 dw 



(18) 



--l/c 



Thus knowing the time dependence of qj and Lij, de- 
termines the time dependence of the coefficients of the 
Joukowski map, and hence the evolution of the droplet. 
From the above equations we obtain 



qi = = 

c 



Mi(0)e- t/T = 



r 

02 = -- 

c 

uq 



U C 



Ui, 



(19) 



Li 2 (0)e- t/T =u 



up 
(c 2 - l) 2 



where qj and fij(0) are determined by the initial shape 
of the droplet, and without loss of generality we choose 
qi = 0. The evolution that follows from the above equa- 
tions, with the initial conditions r(0) = 1/100, Mo(0) = 
— 1/2, ui(0) = —5/12, and c(0) = 6/5, is demonstrated in 
Fig. 3. It describes a droplet pinching of the puddle. The 
above equations hold up to the point where the droplets 
separate. Beyond this point the two droplets continue to 
evolve independent of each other. Namely in this par- 
ticular case the droplets remain circular while their area 
decrease exponentially in time. 

Let us now describe the evolution of droplets whose 
boundaries are described by conformal maps of the 
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formh 



z(w) 



rw 



w 

W ; 



(20) 



with the conditions \u>j\ > 1 and a j = 0- As m the 
previous case, to extract the time dependence of the co- 
efficients, r, a,, and wj we shall identify the singularities 
of the Schwarz function, and require them to agree with 
the time dependence dictated by ([5]) and (fT5|) . 

From (fT4|) and (|20|) it follows that the singular points 
of S(z) in I?, denoted henceforth by qj, are the points 
for which qj = z{\/wj). The behavior of the Schwarz 
function near these points is given by 



S {qj + 5z) = z 



Sz 



^ln((5z), (21) 



The conditions 



and therefore S(z) ~ . otj ln(z — g>j) 

= 0, ensures that no branch cut 



> 1 and £\a 
crosses the droplet boundary, and from (|15p . we obtain 
tk = —J2j a jQj +1 /k, for k > 0. Thus, the uniform ex- 
ponential decay of the harmonic moments implies that 
otj(t) — aj(0)e^*/ T , while qj are constants of motion. 
The equation for to, is obtained by deforming the con- 
tour of the integral (|15p around the pole at infinity. The 
singular behavior of the Schwarz function in this region 
is 



S(z) 



rr \ ~t _ r 
/ a j— — 

Z Z — ' 1/1- 7. 



r 

WiZ 



oo. 



(22) 



and therefore the zeroth harmonic moment is to = l r | 2 — 
r J^n tt- This relation provides the additional equation 

^ — 'J W j 

needed in order to determine the coefficients of the con- 
formal map (|20[) : 



q j = z(wj 1 ), o i (0)e- t / T =a J -(t), 

Otj(t) 



(23) 



h(0)e 



-t/r = M 2_ r 



E 



W j 



Two examples of the evolution generated by confor- 
mal maps of the type (|20|) are depicted in Fig. l(lo|. 
The left-hand panel shows the evolution of a droplet 
contracting to a star-shape form were each one of its 
arms is similar to the Saffman- Taylor finger in a channel 
geometry ll( . In reality these water fingers eventually 
break into small droplets due to surface tension. Nev- 
ertheless, they have been seen in numerical simulations 
of related problems [1 21. and their traces are observed in 
the experimental data[3j. These studies also reveal that 
a rather stable pattern, that commonly appears in these 
systems, is a dry patch finger split in the middle by a nar- 
row water " pipe" , similar to the result we obtained shown 
in the right panel of Fig. 1. The common appearance of 



these fingers, called "doublons", in volatile thin films is 
in a marked difference with the Saffman- Taylor evolution 
in a radial geometry, where tips split repeatedly. 

It should be noted that Eqs. ((3]) and (j4|) also have 
solutions which evolve into cusp singularities as well as 
tip splitting, unobserved in the experimental data. An 
additional information is needed in order to select the rel- 
evant solutions. This information comes from the local 
behavior of the fluid near the droplet boundary. There 
is evidence that Rayleigh instability generates inhomo- 
geneous distribution of the liquid along the rim of the 
droplet 0. It is expected that this behavior selects the 
relevant solutions, such as the finger shown on the right- 
hand panel of Fig. 1. This problem is left for further 
studies. 
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